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Abstract 

We study various self-avoiding walks (SAWs) which are constrained to lie in the upper half¬ 
plane and are subjected to a compressive force. This force is applied to the vertex or vertices 
of the walk located at the maximum distance above the boundary of the half-space. In the case 
of bridges, this is the unique end-point. In the case of SAWs or self-avoiding polygons, this 
corresponds to all vertices of maximal height. We first use the conjectured relation with the 
Schramm-Loewner evolution to predict the form of the partition function including the values 
of the exponents, and then we use series analysis to test these predictions. 
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1 Introduction 

Self-avoiding walks (SAWs) were initially introduced [9, 28] as a model of long linear polymer 
chains. Since that time they have been explored both as a polymer model in statistical mechanics 
and as an independent problem of interest to combinatorialists and theoretical computer scientists. 

In the simplest case, one considers walks on the edges of a lattice, starting from a fixed origin 
and forbidden from visiting a vertex more than once. We will denote by c n the number of such 
objects with exactly n steps. A simple subadditivity argument [13] shows that 

lim — logc n = f3 := inf — logc n . 

n —¥oo ft n—>00 Tb 

The quantity e 13 is known as the connective constant, and depends on the lattice in question. 
The connective constant is unknown for all but one regular lattice, though numerical estimates have 
been computed for a number of cases. Of interest here is the current best (nonrigorous) estimate 
for the connective constant of SAWs on the square lattice [5], 

= 2.63815853035(2). 
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The one exception is the honeycomb lattice, for which it was predicted in [27] and proved in [7] 
that = y/ 2 + y/2. We will set 


and note that ql/ n -A 1. We also set 


Qn — e ^ c n , 


<? 2 n 

Pn — / \ 2 > 

(< 7± 2 * 


( 1 ) 


which is the probability that two SAWs of length n can be concatenated to give a SAW of length 
2 n. We will be primarily interested in the square lattice Z 2 , but we first review what is known 
about Z d . 

In the trivial case d = 1, c n = q n = 2 ,p n = 1/2 for all n > 1. If d > 5, it is known [14] 
that p n —> poo > 0 ; this is suggested by (but harder to prove than) the fact that two independent 
ordinary walks in five or more dimensions have a positive probability of no intersection. Little is 
known rigorously about q n for dimensions 2,3,4. It is widely believed that 


q n ~ 


An Yd 1 
A (logn ) 1 / 4 


d = 2,3 
d = 4. 


Here A is a lattice dependent constant 1 , 72 = 43/32 (this was first predicted by Nienhuis [27] and 
also follows from SLE computations, see Section 4); 73 = 1.156957 ± 0.000009 [3]. For d = 4, see 
[1] for closely related results on a similar model. In particular p n x n ~ 11//32 for d = 2. For ordinary 
walks without self-avoidance, it is known [ 20 , 22 ] that p n x n “ 5 / 8 where p n here is defined to be 
the probability that the set of points visited by two simple random walks of n steps starting at 
adjacent points are disjoint. (In our notation, a n ~ b n means that a n /b n —> 1 and a n x b n means 
that there exists c such that a n /c < b n < ca n for all n sufficiently large.) 

The average size of an n-step SAW is often measured by the mean squared end-to-end distance 
IE n [|cj n | 2 ], where E n denotes expectation with respect to the uniform probability measure on walks 
of length n starting at the origin. It is widely believed that 


E n [|Wn| 2 ] ~ 


A n 2 ", d^ 4 

An (logn) 1 / 4 , d = 4 


where v is a dimension dependent exponent. In the trivial case d = 1, we have n = 1 and for 
d > 5 this has been proved [14] with v = 1/2. In two dimensions it is predicted (see [27] and 
the SLE derivation below) that v = 3/4, while in three dimensions the best current estimate is 
v = 0.587597 ±0.000007 [4], 

For the remainder of this paper, we set d = 2 and write Z 2 = ZxiZ. The theoretical (although at 
this point not mathematically rigorous) understanding of two-dimensional SAWs has been deepened 
by considering the implications of assuming that the walks have a conformally invariant scaling 
limit. Rather than working with walks of a fixed number of steps, it is more useful to consider the 
measure on SAWs of arbitrary length that gives measure e~^ n to each SAW of n steps, with any 
starting vertex. We will write Q for this measure, that is, for any set V of SAWs, 


OO 

Q(V) = ^e-P n c n (V), 

n =0 

1 Thoughout this paper we will use A for a lattice dependent constant; however, its value will change from line to 

line. However, when we write Ai, A2 ,..., these values will stay fixed. In all cases, the values of these “nonuniversal” 

constants are not important to us. 
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where c n (V) is the number of walks of length n in the set V. In particular, if V n is the set of SAWs 
of length n starting at the origin, Q(V n ) = q n . 

For any lattice spacing 1 /IV, we can write Qn for the measure Q viewed as a measure on 
scaled SAWs on the lattice IV -1 Z 2 . Then (see Section 4 for more detail) if it is assumed that 
the measures Qn have a limiting measure that is conformally invariant, then this limit must be 
a version of Schramm-Loewner evolution with parameter 8/3 which we denote by SLE 8/ / 3 . (There 
is a one-parameter family of SLE curves, but there is only one value 8/3 that gives a measure on 
simple curves satisfying a “restriction property” which would have to hold for a scaling limit of 
SAWs. For this reason we will only discuss SLE 8 / 3 in this paper. The values for exponents below 
such as b,b are particular to the parameter 8/3.) The analogues of the exponents v and 7 can 
be computed for SLE 8 / 3 and give v = 3 / 4,7 = 43/32. These are mathematically rigorous results 
about SLE 8 / 3 although it is still only a conjecture (with strong theoretical and numerical evidence) 
that it is the limit of the measure Q. 

In this article we consider the enumeration of several different types of SAWs constrained to 
lie in the upper half-space of the square lattice, weighted with a Boltzmann weight corresponding 
to the maximum distance above the boundary of the half-space reached by the walks. When this 
weight is less than one, walks which step far away from the boundary are penalised, and one 
can thus view this system as a model of polymers at an impenetrable surface, subject to a force 
compressing them against the surface. 

2 The model 

Let W denote the set of SAWs in Z 2 = Z x iL starting at the origin and let W + be the set of 
half-plane walks in W, that is, walks that stay in 

H := {x + iy £ Z 2 : y > 0}. 

We let W n , W+ denote the set of such walks of length n. 

If oj = [0, uj\ ,..., oj n ] £ W + we write 

h{oj) = max{Im(uj) : j = 0,1,..., n}, y(w) = Im(o; n ). 

We call a walk oj = [wq, •.. ,w n ] £ W+ a bridge if h(pj) = y(oj). The utility of bridges lies in the 
fact that they can be freely concatenated to form larger bridges (with the addition of an extra step 
between). We write B n for the set of bridges of length n and B for the set of all bridges. We write 
(•}, (•)+, (-) b , (•)„, (•)+, (•}(( for integrals (in fact, sums) with respect to the measure Q restricted to 
the sets W, W + , B, W n , W^,B n , respectively. For example, 

(l)n = ^2 1 ' < 3( w ) = Cne ^ = V™- 

Note that (1) = (1) + = (l ) 6 = 00 , as the corresponding generating functions diverge at j3. 

A walk oj = [cuo ..., o: n - 1 ] £ W n _i is called a polygon if \oi n -\ — wo| = 1, that is, if the first and 
last vertices of oj are adjacent. An edge can be added between ojq and cu n _i to form a simple closed 
loop. We say that oj has length n. Let V be the set of all polygons and V n the set of polygons 
of length 71 , with V + and the analogous sets restricted to those polygons staying in HI. Then 
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let (-) p , (•)«, (-) p+ , (")n + be the integrals with respect to the measure Q restricted to V, V n . V + , V+ 
respectively. 

A polymer model which has been considered in the past [11, 15, 19, 30] is that of polymers 
terminally attached to an impenetrable surface, with an external agent exerting a force on the non- 
attached end of the polymer, in a direction perpendicular to the surface. This reflects real-world 
experiments where polymers are pulled using optical tweezers or atomic force microscopy [31]. This 
can be modelled using half-plane SAWs with a Boltzmann weight associated with the height of the 
end of the walks above the surface. The partition function of such a model is 

CT\u) = <e-^)>+. 

One can interpret u as the reduced pulling force; when u < 0 the force is pulling up away from the 
surface, and when u > 0 the force is pulling down towards the surface. As we can see below, it is 
predicted that for fixed y, as n —> oo, 

= »- 19 / 1 V 5/28 , 

and hence the integral is dominated by the small values of y which gives the following. 
Prediction 1. For each u > 0, there exists a constant A = A[u ) such that 

An ~ 19/16 . 

It has recently been shown [2] that this model undergoes a phase transition at u = 0; that is, 
the free energy 

k(u) := lim - log C and (u) + 0 

n—> oo Ti 

is non-analytic at u = 0, being equal to fi for u > 0 but strictly greater than j3 for u < 0. This is 
the transition between the free and ballistic phases of self-avoiding walks. 

We will consider three similar quantities 

CT» = (e- uHu) )^, (2) 

Bn{u) = (e~ uh ^) b n = {e~ uy ^)l (3) 

P n ma » = {e~ uh ^y n + (4) 

in the u > 0 regime. We will use SLE to give predictions for the asymptotic values and then will 
analyse exact enumerations to test the predictions. (For polygons, a numerical analysis has been 
conducted elsewhere [12]). In this case the integral does not concentrate on the lowest order terms, 
and there will be a stretched exponential decay. We state the prediction now. We note that the 
constants Ai and A 2 do not depend on u. 

Prediction 2. There exists Ai, A 2 > 0 such that for each u > 0, there exist u-dependent constants 
A+, A b u , A y+ such that 

C™ x (u) ~ A+ n 3/16 exp{-Aiu 4/7 n 3/7 }, 

B n (u) ~ A b u n~ 13 / 28 exp{—A lU A/7 n 3/7 }, 

P” ax (u) ~ A p u + n~ n/ 7 exp{-A 2 u 4/7 n 3 / 7 }. 
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3 A simpler model: ordinary random walks 


Here we consider an analogous problem for simple random walks that can be solved rigorously. Let 
(■)n ,simp denote expectations with respect to simple random walks of length n started at the origin 
restricted to stay in the upper half-plane HI. In other words, 

A)n’ sim P = 4 _n Y (v) 
v 


where the sum is over all nearest neighbour (not necessarily self-avoiding) paths of n steps starting 
at the origin with 77 C El. We write h(r]) for the maximal imaginary component of ij. We will give 
the asymptotics of ( e - M M w )}+> smip for u > 0 as n —> 00 . 

Let c(n, h: r, y) denote the number of simple random walk paths of n steps starting at height 
r; ending at height y; whose height stays between 0 and h — 2 for all times. Then by viewing the 
imaginary part as a one-dimensional random walk killed when it leaves the interval [0, h — 2], we 
can see that 4~ n c(n, h; r , y ) = y) where { Jh{r, y) : 0 < r,y < h — 2} is the symmetric matrix 
given by 

J h(r,r) = ~, J h (r,r± 1) =-. 

The eigenvalues and eigenfunctions of Jh are well known, and can be computed using a Fourier 
series (sum) over {0,1,... ,h — 2}, and simple trigonometric identities. If Vj denotes the vector 


with components 


sm 


j(k+l)ir 

h 


J k=0,...,h-2 


, then 


1 1 


_ _ j 7f 

Jh Vj = Aj Vj, where A j = ^ ^ cos ( — 


By diagonalizing J, we can see that 


h — 1 


4”(n v) = T X 


3 =1 


1 1 


- + - COS 1 — 

h 


jit 


sm 


(r + 1)jtt 

h 


sm 


(y + l)j7r 

h 


If h, n —> 00 with h 2 <C n, then the j = 1 term dominates the asymptotics and 

Jh(r,y) 


2 

h 

2 

h 


1 1 

—I— cos 
2 2 


sm 


(r + l)7r 
h 


sm 


(V + 1)tt' 
h 


exp 


In particular, 


f mr 2 1 

sin 

(r + l)7r 

sin 

~{y + i)tt" 

\ 4/7 2 / 

h 

L 

h J 

2n 

J U7T 2 I . 

"(y + 

1)7T 


" 7? exp 

\ 4/i 2 J Sm 

h 




Let F n (h) denote the probability that a simple random walk starting at the origin in Z 2 up to 
time n stays in HI and has maximum imaginary component less than or equal to h — 2. Then, if 
1 </l< y/n, 


h -2 

F n (h) = J2 J h( 0,2/) 

y =0 


h 


exp 


TITT 2 1 

4P"J ' 


( 5 ) 
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Note that F n (h) — F n (h — 1) gives the probability that the maximal imaginary component is exactly 
h— 2. While we used the exact solution to determine these asymptotics, one can give a short heuristic 
to explain why the answer should be of the form ch _1 exp{— u(n/h 2 )} for some c, u. In the first h 2 
steps, the probability that a walk stays in the upper half plane is 0(/i -1 ) by the “gambler’s ruin” 
estimate for one-dimensional walks. After this, every time the walk moves h 2 steps, the chance 
that its imaginary component leaves [0, h — 2] is strictly between 0 and 1. If we call this probability 
1 — e~ a , then the probability that the walk stays in for n steps should look like [e~ a ] n ^ h2 . As 
h —> oo, the exponential constant a can be computed as an eigenvalue for Brownian motion. 

We see that 


(e ~ uh )+> simp = ^ e~ uh [F n (h + 2) - F n (h + 1)] 

h =0 

oo 

= e 2u Y J Z~ uh [F n (h)-F n (h- 1)] 

h =2 

oo 

= e 2u {l - e- u )^2e~ uh F n (h) 
h= 2 

oo ( 2 ^ 

~4 e 2 «(l-e-«)£frV“‘exp{-^U, 

h =2 J 

where we have used (5) for h <C \fn. For larger h, both e~ uh F n {h ) and h~ 1 e~ uh exp{— ra7r 2 /4/r 2 } 
decay quickly to 0. Then using (18) in the Appendix, 



where A u , X u can be derived from (18): 

A u = 2 8 / 3 • 3" 1 / 2 • vr 1 / 6 • u~ 1/3 • e u (e u - 1), 
X u = 2 -4 / 3 • 3 • 7r 2 / 3 • u 2 / 3 . 


4 Scaling limit 

The structure of this section is as follows. In Section 4.1 we briefly review Schramm-Loewner 
evolution and demonstrate that, under the assumption that SLE 8 / 3 is the scaling limit of two- 
dimensional SAWs, it can be used to predict the asymptotic behaviour of the Q-measures of SAWs 
under various types of restrictions. In Section 4.2 we turn to the specific problem of compressed 
SAWs, which can be phrased in terms of walks restricted to a horizontal strip of fixed height, and 
together with Lemma 1 in the Appendix derive the first two results in Prediction 2. In Section 4.3 
we consider polygons restricted to a horizontal strip. These require a slightly different approach, 
and we use results derived from restriction measures and large deviation theory to obtain the third 
part of Prediction 2. 
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4.1 Schramm-Loewner evolution 

Here we review predictions for the measure Q using the Schramm-Loewner evolution (SLE) in [24] 
and extend them to incorporate information about the natural parametrization [25, 26]. We start 
by explaining precisely what is meant by the scaling limit of SAWs. The starting point is to assume 
the existence of the mean-square displacement exponent v and two scaling exponents b and b that 
we discuss below. No a priori assumptions about the values of v, b, b are made. It is convenient 
to view a SAW uo = [coo? • • • ,Wn] as a continuous time process uj : [0, n] — > Z 2 with uo(k) = and 
otherwise defined by linear interpolation. The exponent v is expressed by stating that the fractal 
dimension of the paths in the scaling limit is 1/v. For every lattice spacing 1/N and c > 0, and 
every SAW lo = [wo, wi, • • •, uj n ], we define the scaled walk oj n on the lattice iV -1 Z 2 by 

uj N (t) = N- 1 u(tN v /c), 0 <t< cn/N u . (7) 

The measure Q induces a measure on scaled paths in N~ l Z 2 . Note that we do not change the 
measure except that we view it as a measure on scaled paths. If uj has length n, then the time 
duration of the curve uj n is cn/N v . We discuss below what value to choose for c. 

Let HI denote the upper half-plane as before. If z E HI and n is an integer, let Vo }nz denote 
the set of SAWs in Z 2 starting at the origin and ending at nz. (We are assuming for notational 
ease that nz is a lattice point; if not, choose the closest such point.) Let V^ nz denote the set of 
half-plane SAWs in Vo,nz- The various types of restricted SAWs considered in this section can all 
be regarded in terms of Vo,nz and V^ nz , and so the limiting behaviour of the Q-measures of these 
sets is central to our predictions. 

• Assumption: There exists a boundary scaling exponent b, an interior scaling exponent b, 
and continuous functions iJj, on HI such that as n —> oo 

Q[Vo,m] ~ ifj(z)n~ 2b . (8) 

Ql v o%] ~ ^+(A) n~ b ™~ b - (9) 

The scaling exponents 6, b , come from the power law behavior in the total mass (partition function) 
of the measure Q. 

The scaling limit of the SAW is a collection of measures {/u(0, z), //h(0, z)} that are the limit 
of the measures n 2b Q and n b+b Q restricted to the sets Vo tnz , I7 0 + n2 , respectively. These are finite 
measures on simple curves from 0 to z in C or HI, respectively. They have total mass (normal¬ 
ized partition function) i/j(z),i/j + (z) and one can get probability measures {/x # (0, z), //^(0, z)} by 
normalization. 

More generally, take z, w to lie in the interior of HI (i.e. not on the boundary of the half-plane), 
and define V nz/nw to be the set of SAWs in Z 2 starting at nz and ending at nw. Then the measure 
fj(z, w) is the limit of the measure n 2b Q restricted to the set V nz . nw . Then if D is an open subset of 
C, and w are distinct interior points, we can define w ) to be the measure /jl(z, w) restricted 

to curves that stay in D. Equivalently, it is the scaling limit as above of SAWs from nz to nw that 
stay in nD. This collection of measures satisfies the restriction property: if D C D' and z, w E D, 
then ii£)(z,w) is /j,£>'(z,w) restricted to curves that stay in D. Similarly, if z E dD and w E D we 
can define the measure /id(z,w) provided that the boundary is sufficiently smooth at z (there are 
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lattice issues involved if the boundary of D is not parallel with a coordinate axis (see [18]), but we 
will not worry about this here.) 

In addition to the restriction property, we will assume that the measures hd(z,w) satisfy two 
further criteria. First, we assume that the limit is in some sense conformally invariant: 

LL 

• The probability measures z , w) considered as a measure on curves modulo reparametriza- 

tion 2 are conformally invariant. 

-LL 

Given this, one can define the probability measures /i/)(z, w) even for boundary points at which 
the boundary is not smooth. This allows one to write down another property that the scaling limit 
of SAW should have. It can be considered a property of curves up to reparametrization. 

• Domain Markov property. Suppose an initial segment w[0,f] is observed of a curve from 
/j^(z,w). Then the distribution of the remainder of the curve is that of /r^,(u;(i), z) where 
D' is the slit domain D \ w[0, t]. 

The three properties: conformal invariance, domain Markov property, and restriction property 
characterize the scaling limit in that there is only one family of measures on simple curves that 
satisify all three. The measure is called SLE 8 / 3 , the Schramm-Loewner evolution with parameter 
8/3. (There are SLE measures with other parameters but they do not satisfy the restriction property 
which a scaling limit of SAW would necessarily satisfy.) 

The critical exponents for SAW can be deduced theoretically (but not at the moment math¬ 
ematically rigorously) from mathematically rigorous theorems about SLE 8/ / 3 . First, the fractal 
dimension of SLE 8/ / 3 paths is 4/3. This was first proved as a statement about the Hausdorff dimen¬ 
sion, but for our purposes, it is more useful to think of it in terms of the A/S-Minkowski content. 
If X is a compact subset of C, then the 4/3-dimensional Minkowski content of X is defined by 

Cork (X ) = Cont 4 / 3 (A) = lime -2 / 3 Areajz : dist(z, X) < e}, 

provided that the limit exists. If co : [0,T] — > C is a curve, then we can view Cont(u;[0, t]) as the 
“4/3-dimensional length” of oo[0,t]. It has been shown that for the SLE 8 / 3 measures, the function 
t i—)• Cont(cu[0, t]) is continuous and strictly increasing. For this reason we can parametrize our 
SLE 8 / 3 paths so that at each time t, Cont(u;[0, t]) = t. This is called the natural parametrization 
and is the 4/3-dimensional analogue of parametrization by arc length. For the remainder we will 
assume that we have parametrized our curves in this way. We conjecture that we can choose c in 
(7) so that the curves in the scaling limit have the natural parametrization, and for the remainder 
we assume we have chosen this c. 

If uo : [0,T W ] —> C is a curve taking values in a domain D and /:£)—>• f(D ) is a conformal 
transformation, then we define the image curve / o oo to be the image with the parametrization 
adjusted appropriately. To be precise, the time to traverse /(cj[0,f]) is 

/Vms»i 4/3 *- 

Jo 

2 That is, if two curves 7 : [0,1] —> C and 7 : [0,T] — > C trace out the same path (but over different times), they 

are viewed as the same curve. 



For example, if f{z) = rz with r > 0, then the total time to traverse the curve is multiplied by 
?’ 4 / 3 ; this is the scaling property of a “4/3-dimensional length”. If /x is a measure on curves in D, 
then we define the measure / o ^ by 

f o ii(V) = n{u : f ou £ V}. 

We can now define the SLE 8 / 3 measures xu)} as a family of measure with the following 

properties. Here we can choose z,w £ D (whole-plane SLE 8 / 3 ); z £ dD,w £ D (radial SLE 8 / 3 ) or 
z,w £ (9.D (chordal SLE 8 / 3 ). In the case of boundary points, we assume that dD is locally analytic 
around the points. The conformal covariance rule is 

fon D (z,w) = \ f'(z)\ b * \f\z)\ bw /x/(u)(/(«),/(«;)), 

where b z and b w take values b or 6 , depending on whether z and w are boundary points or interior 
points. We can write 

Vd(z, w) = $d{z, w ) n*(z, w ), 

where $>d(z,w) is the total mass (normalized partition function) of the measure hd{z,w). The 
covariance rule can then be stated as conformal invariance of the probability measures and a scaling 
rule for the partition functions, 

$d(z,w) = \ f\z)\ bz \f{z)\ bw $f {D )(f(z)J(w)). 

Rigorously, we can show that 



and the chordal, radial, and whole plane partition functions are each determined uniquely up to 
a mulitplicative constant by the scaling rule and the restriction property. For one choice of these 
constants, we conjecture that we get the scaling limits of SAW as above. We will assume this value 
of the constant (we do not know, even nonrigorously, this value). 

In order to relate the values of SLE 8 / 3 to SAWs of a given number of steps, we must consider 
the time duration T u of the paths. While it is known that the time duration of SLE 8 / 3 paths can 
be given by the Minkowski content, the next statement is currently only a conjecture about SLE 8 / 3 . 

• The measures ^d{z,w) can be written as 

Hd{z, w, t) dt , 

where for each t, no{z,w,t) is a (strictly positive) finite measure on curves u from z to w 
with TL, = t. We can write 

jj 

H D (z,w,t) = $d(z, w,t) n]j(z,w;t), 

n 

for a probability measure /j^(z,w;t) with $£>(z,w,t) continuous in t. 

In other words, with respect to the probability measure fi^(z,w), the random variable T u has a 
strictly positive density given by &d(z, w, t)/$D(z, w). 

We now use this to compute (nonrigorously) the exponent 7 for SAW in Z 2 . For the remainder 
of the paper, if we establish that two sequences a n and b n satisfy a n x b n , we will often conclude by 
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conjecturing that there exists c such that a n ~ cb n . We will not justify this last step, and, indeed, 
in most cases we are just assuming that things are “nice”. 

Let m n denote the Q measure of the set of SAWs starting at the origin of length between n 4 / 3 
and (2n) 4//3 . Then m n x n 4//3 q n 4/3. Typically, these walks should go distance of order n and hence, 
we would expect that the Q measure of the set of SAWs starting at the origin of length between 
n 4 / 3 and (2n) 4 / 3 whose endpoint is distance between n and 2n is also comparable to m n . There are 
of order n 2 possibilities for this endpoint, and the number for each of the possible endpoints should 
be comparable, so we would expect that the number of SAWs starting at 0 ending at n of between 
n 4 / 3 and (2n) 4 /' 3 steps is comparable to m n /n 2 . By the scaling assumption, this is comparable 
to n~ 2b = n -5 / 24 . Hence m n x n 2 n~ 5 / 24 = n 43 / 24 = [n 4 / 3 ] [?r 4//3 ] n / 32 . This gives the prediction 


q n 4/3 X [„</3] 11/32, 


Prediction 3. There exists A such that 


(l)n = Qn ~ An 11/32 . 

The argument to compute m +, the Q measure of SAWs starting at the origin of length n 4 / 3 to 
2 n 4 / 3 restricted to the half-plane, is similar. The only difference is that we now have one boundary 
point. Hence, instead of n~ 2b we have n~^ b+b ^ = ?r -35 / 48 , and x n 2 n -35 / 48 = n 61//48 = 
[n 4 / 3 ] [n 4 / 3 ] -3 / 64 . 


Prediction 4. There exists A + such that 


(1)+~ A+n" 3 / 64 . 

Finally, to compute m^, the Q measure of bridges starting at the origin of length n 4//3 to 
2 n 4 / 3 restricted to the half-plane, we now have two boundary points. Hence, m^ x n 2 n -5 / 4 = 
[n 4 / 3 ] [n 4 / 3 ]" 7 / 16 . 

Prediction 5. There exists A b such that 


(1 }^A fc n- 7 / 16 . 


Predictions 3 and 4 imply that the probability that a SAW of length n is a half-plane walk is 
comparable to n _25,/64 , and Predictions 4 and 5 imply that the probability that a half-plane SAW of 
length n is a bridge is comparable to n -25//64 . (It is not surprising that these are the same - in the 
first case, we take SAW and insist that no vertex lie below the starting point; in the second case, we 
take a half-plane SAW and insist that no vertex lie above its final point. Symmetry suggests that 
the fraction of walks we keep ought to be roughly the same in both cases.) Hence the probability 
that a SAW of length n is a bridge is comparable to n -25 / 32 . 

We finish by mentioning two other well-known conjectures. Let be the indicator function 
that the endpoints of u are distance one apart. Then observe that (l)^ = ( $ ) 2 n-i and {l) p 2 l = 

Prediction 6. There exist A, A' such that 

<*> 2 n-l = (l>2n ~ An~ 312 , (<&>+_, = <1}£ ~ A' n^ 2 . 
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For the first relation, we note that the measure of SAWs of 2 n steps whose middle vertex is 
the origin is comparable to n 11 / 32 . The endpoints of these walks are distance n 3 / 4 away and hence 
there are of order n 3 / 2 possibilities for each of the endpoints. We therefore get a factor of n ~ 3 / 2 
which represents the probability that the endpoints of the 2n-step walk agree. We also require the 
absence of any intersection between the first and second halves of the walk near the initial/terminal 
point. This gives another factor that is comparable to ?r~ 11//32 . Therefore, the measure should be 
comparable to 

n 11 / 32 n -3 / 2 n -11 / 32 = n -3 / 2 . ( 10 ) 

One may note that the exponent 7 cancels in this calculation — one did not need to know its 
value. This shows that SAPs can be easier to analyze than SAWs since there are no endpoints. We 
will use the basic principle of (10) several times below, so we state it. This idea extends to other 
dimensions, so we state it a little more generally. 

• The probability that the endpoints of a SAW of length 2n — 1 are distance one apart is 
comparable to p 2 n n~ du , where p 2 n is as defined in (1). In particular, if d = 2 it is comparable 

to 77. 14 / 32 77, 3 / 2 = 77, 59/ 32 ^ 

For the second relation in the prediction, we start with a SAP and consider the various loops 
one gets by translating the root around. Typically there are only 0(1) such vertices with minimum 
imaginary part, and hence the probability that a SAP of length 2 n is a half-plane SAP is comparable 

to 77 , 1 . 

Now let q+ h denote the measure of the set of SAWs of length n that stay in the upper half-plane 
and that start at ih. If h = 0, then q+ 0 = qf and if h > n 3 / 4 , then q+ h x q n . We will consider 
the case 0 < h <C n 3 / 4 , and define K by n = I\ / 1 4 / 3 . Then a SAW of n steps can be viewed on 
a “mesoscopic” scale as a SAW of K steps where each of the steps is a little SAW of / 1 4 / 3 steps. 
Using this perspective, we see that we would conjecture that the probability that a SAW of Kh 4 / 3 
steps starting at ih stays in the half-plane is comparable to AT -25 / 64 . Hence we conjecture that qif, 
is comparable to q n AT -25 / 64 which is comparable to qif / 1 25 / 48 . To phrase this in our notation, let 
Sh(u}) be the indicator function that Im(cjj) > —h for all j. 

Prediction 7. There exists A such that if 1 <C h <C n 3//4 , then 

(Sh)n ~ An -3 / 64 ft 25 / 48 . 

If z E H, let q n (z),qf(z) denote the measure of SAWs and half-plane SAWS, respectively, of 
length n starting at the origin ending at z. If z = x + iy with x x n 3//4 , y x n 3//4 , then we expect 
that q n (z) x n -3 / 2 q n , qif(z) x n -3 / 2 qd. Now suppose x x n 3//4 and 1 <C y <C n 3 / 4 . Then we view 
a SAW from 0 to y as the concatentation of two walks of length n /2 — a half-plane walk starting at 
the origin and (the reversal of) a SAW starting at x + iy that stays in the half-plane. We therefore 
guess that the measure should be comparable to n ~ 3//64 [n _3//64 y 25 / 48 ] n ~ 59 / 32 = n -3i/i6 ^25/48_ jj ere 
the 59 / 32 = n ~ 3 / 2 n , -11 / 32 as in the principle above. For a fixed y, there are of order n 3 / 4 possible 
values of x. By summing over these we get the following prediction. 

Prediction 8. There exists A such that if 1 <C h <C n 3//4 , then 

(1 [y{u) , h) r~An-'0^h^. 
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4.2 SAWs restricted to a strip 

The predictions in the last subsubsection were made assuming that the the maximal imaginary 
displacement of the bridge or half-plane walk was typical for the number of steps. We now consider 
the case where this maximal height is much smaller. We will fix a height h and consider SAWs 
restricted to the infinite strips 

Sh = {x + iy E Z x iZ : 0 < y < h}, 

„ h f „ „ h h) 

S h = S h - - i= lx + iy £ZxiZ ■. -~<y < - \ . 

We will write Sh,S^ for the indicator functions that oj E Sh,S respectively. We will consider 
SAWs in Sh of length n = K /i 4 / 3 where 1. We write any such SAW w as a concatenation of 

K SAWs of length h 4 / 3 : 

OJ = OJ 1 © • • • © 0JK- (11) 

Let us choose a lattice space of 1 /h and write the scaled walks as 

00^ = UJ]‘ © • • • © CJjk. 

These walks live in the strip 

D = {x + iy E C : 0 < y < 1}. 

The walk oj h can be viewed as a polymer consisting of K monomers cjj each of length of order 1. 
On this scale w is a one-dimensional (in some sense, weakly) self-avoiding walk of K steps. Hence, 
we would expect the measure of such walks to decay like e~^ lK for some /3±. Moreover, we expect 
no smaller order corrections in K (that is, the measure is asymptotic to A\ e~^ lK for some A\ ) 
and for the typical walk to look like a straight line of length cK. This is the basis for the following 
predictions. We will assume that h —> oo and Kh —> oo and write n = Kh h 4 / 3 . 

Prediction 9. There exist f3\ and A±, A 2 , A 3 , A 4 such that the following holds. 

(Sh) n ~ Ai exp |-^in/i _4/3 | h 11/24 , 

(S h )t ~ A 2 exp |-/3in h _4/3 | h _1/16 , 

(1 {h(u)=h})Z = {Sh)n - (Sh-i)n ~ ^3 exp|-/3m/i“ 4/3 } [n h~ 7/3 ] h~ 1/w , 

{^-{y(u)=h)}) b n ~ A 4 exp|-/3m/l _4/3 } h- 19 / 12 . 

The first two relations follow from Predictions 3 and 4 respectively, where the factor of e~^ lK 
accounts for the height restriction and we have used n = Kh 4 / 3 . The third is merely the L-derivative 
of the second. The fourth follows from Prediction 5. We note that /i~ 19 / 12 = /i -7 / 12 hr 1 . The extra 
factor of comes from the fact that we are specifying the exact height. Using (18), we see that 
as n —> 00 , 



Au~ 5 l 7 n 3 / 14 exp j — Aj. u 4 ^ 7 n 3 / 7 } , 
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rsj 



where 


Ai = 7 • 3 -3 / 7 • 4 _4//7 ■ P^ 7 , 

and A does not depend on n or u. Similarly, using (18), we see that there exists ci, c 2 (independent 
of u) such that with the same value of Ai, 


{e-”' ,( ")>+ = £e-”*(l w „,, w >+ 

h= 1 
00 

~ ex p{-ftn ; r 4/3 } [nr 7 ^- 1 / 16 


h =1 


f 

Jo 


A%n / x 


,-115/48 


exp 


Pi n 
4/3 


X' 


+ ux > dx 


ci n 3/16 u 5/16 exp{-A!U 4 / 7 n 3 / 7 }, 


( e -uh ( W ))6 


h=l 


£e -“^ 4 exp {-/3m h- 4 / 3 } h - 19 / 12 
h=l 


= A4 1 X 19/12 exp{-(^+ ux) 
= C 2 n _ 13 // 28 u -1 / 28 exp{—A^/V/ 7 } 
We have thus obtained the results of Prediction 2. 


dx 


4.3 Polygons restricted to a strip 

We finish this section with some predictions for self-avoiding polygons in a strip. 

Prediction 10. There exist ft 2 o,nd A$,Aq such that the following holds. 

{S* h $)n-1 = (StZ ~ A 5 exp {-^nh" 4 / 3 } / 1 ” 17 / 6 , 

= (S'h)P ~ A 6 exp {-ftnr 4 / 3 } hr 25/b . 

We will take a slightly different approach here. For any 0 < y,y < h and integer N > 0, let 
JCh(y, y, N] r) denote the set of SAPs u = [wq, ..., uj 2 n -i,^ 2 n = wq] with the following properties: 


UJ 0 = iy , W2n-1 = i{y - 1 ) 


u r = N + iy, u r+ i = N + i(y- 1),, 

1 < Re[wj] < N — 1, j 7 ^ 0, r, r + 1, 2n — 1, 2n. 

w C 5^, 
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and let 


£h(y, y,N) = \J K(y, y, N; r), K hjN = |J K(y, y, N). 
r y,y 

In other words, K-hN is the set of polygons in S/ whose minimal real value is 0; maximal real value is 
N ; that have only one bond on each of the vertical lines Io := {Refz) = 0} and In ■= {Re(z) = IV}; 
and are translated and oriented so that they start at the “higher” point on Jo and that the final 
bond is the unique bond on Iq. Let lCh,N denote the set of polygons in 5/, whose minimal real 
value is 0; maximal real value is IV; and such that the initial vertex is on Iq. Let Jh,Ni Jh,N denote 
the corresponding indicator functions. In the set lCh,Ni the SAP can have multiple bonds on the 
extremal vertical lines Iq,In', however, we predict that typically there are only a few such bonds. 
In particular, we predict that there exists A such that as h, N — >■ oo with N/h —> oo, 

{S* h $j htN )~A(S* h *J htN ). 

If oi e K,h(y,y,N), then we can view oj as two SAWs: oi\ from iy to N + iy and 012 from i(y — 1 ) 
to IV + i(y — 1). The polygon oj is obtained by concatenating oi\ and (the reversal of) 012 adding 
the two extra bonds to make this a polygon. Any oj\,o)2 can be chosen provided that oi\ n 012 = 0 . 
So we need to find the Q measure of the set of pairs of such walks (011,012) with ui\ n 012 = 0 . The 
limiting measure is predicted to be the restriction measure with exponent 2. To be precise it is 
predicted [21, 23] that the measure is asymptotic to 

Ah~ 4 H D (i(y/h),K + i(y/h)) 2 

where Hp denotes the boundary Poisson kernel (derivative of the Green’s function), K = N/h and 

D = Dk = {x + iy : 0 < x < K, 0 < y < 1}. 

The factor /i~ 4 should be viewed as a factor of h~ 2 for each boundary end of the SAP where 2 is 
the exponent of the restriction measure. Using a conformal transformation we see that 

Ho(i(y/h),K + i(y/h )) ~ c sin(y7r//i) sin(y7r/h) e~ Kn . 

If we sum over the number of possible values of (y, y), we see that we get the prediction 

Q(IC h , N ) ~ Q(JC h , N ) ~ h~ 2 e~ 2K \ 

Now consider the set tCh,N of self-avoiding polygons in starting at the origin such that the 
real displacement equals N= K h. To each such polygon, we can find a polygon that visits the same 
points in the same order that starts at a point of minimal real part. Conversely, if we have a walk 
in K-hN and choose a vertex in the middle, we can translate the walk so that vertex is the starting 
point. Since the typical such polygon has length comparable to K Ii 4//3 , there are 0(A'h 4 / 3 ) choices 
for the vertex. We therefore get the conjecture 

Q(t KN ) ~ A I< h i/3 h~ 2 e~ 2Kn = A I< h~ 2/ 3 e~ 2Kn . 

This gives a prediction for a fixed real displacement, but we must convert it to a fixed number 
of steps. Consider a walk in fCh,N chosen from the probability measure given by Q restricted to 
lCh,Ni normalized to be a probability measure. Let T denote the number of steps of such a walk. T 
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is the sum of K random variables, each representing the number of steps in one of the K squares 
of side length h. These random variables are roughly identically distributed and have short range 
correlations, so as an approximation we view T as having the behavior of 

h 4 / 3 [*! + •••+ X K \, 


where Xj are independent, identically distributed. We write a for the mean. Hence, we predict 
that the expectation of T is ah^^K and the standard deviation of T is comparable to /i 4 / 3 K 1 ! 2 . 
Even more precisely, standard arguments from large deviation theory lead us to expect that (at 
least for small r) there is a rate function p(r) with p(a ) = 0 such that the probability that T = 2n 
where n = /i 4 / 3 r K = /i 1 / 3 r N is comparable to 


h^K 1 ' 2 


exp{— p(r) K}. 


Hence the Q measure of polygons in °f 2n steps is comparable to 




/r 4 / 3 A' 1 / 2 exp{ -p(r) I\} K h~ 2 N e~ 2nK x N 1 / 2 h~ 5/2 


exp 


and 


<s; 4>>„ x h 5/2 g N 1 ' 2 exp {- [p 


{ _ K»w*) + 


2vr 


N h 




+ 27T 


Nh~ 


As before, to analyze this sum we start by finding the N that minimizes 


f(N) = 


n 


Since 


f'(N) = p 


n 


h 1 / 3 N 


P \h 1 / 3 N 

+ 2ir — 


+ 27T 


N = [p{r) + 27t] N. 


n 


we find ?’q with 


rop’{r 0 ) - p(r 0 ) = 2vr. 


h ^N f,, (mn) =p{r) + 2r ~ rAr) ' 


We assume that p is smooth and we note that it has a global minimum at r = a where p(a) = 
p'(a) = 0. Therefore, we expect ro > a, and if we set A = p(ro), 

p(r 0 + e) = A + p'(r 0 ) e + 0(e 2 ) = A + 2i: + X e + 0(e 2 ). 

ro 

The terms in the sum are maximized at N n where N n h^ :i r\) = n. Moreover, for k <C N n , 

n 


f(N n + k) = 


P 


h V3 (N n + k) 
k 


+ 27T 


k 2 


(N n + k) 


P r 0 - r 0 — + O —x + 2vr 


A- 


Nn 

(A + 2n)k 

iVn 


+ o 


( k 


V N 2 


+ 27T 


(N n + k) 
(N n + k) 


= /W,) + o(0 
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The key fact about this calculation is that f(N n + k) is comparable to f(N n ) for k 2 < N n , that is, 
for 0 (Nn 2 ) values of k and hence 


N 1/2 

N =1 



and 


+ 2 n 


Nh-^^Nn exp |— /3 2 n/r _4/3 } x fn/h 1 ' 3 ) exp { — /5 2 r?A“ 4/3 } , 


(SI $) n ~nh 17 / 6 exp | -forth 4 / 3 j 


for some /3 2 that we cannot determine explicitly. 

For the second part of Prediction 10, we start by considering a polygon with S? <F = 1 as two 
parts, w 1 and oj 2 . There is a “lower” SAW w 1 of 2/i 4 / 3 steps whose middle vertex is the origin, and 
then cj 2 is an “upper” n — 2h 4 / 3 step SAW from one endpoint of ca 1 to the other that does not 
intersect a; 1 . Using the reasoning for polygons from before, we predict that the probability that w 1 
lies entirely in the upper half-plane is comparable to /i“ 4 / 3 . Hence, we would expect 


(S h $}„ x h~ 4 ' 3 (S* h <h) n 


Finally, we can associate the weight e~ uh with polygons of height h, as we did previously with 
walks and bridges. To do so we must first examine polygons of height exactly h, rather than height 
at most h. 


Prediction 11. There exist constants A 7 and /3 2 (taking the same value as in Prediction 10) such 
that 

(S h <!> l {hH=h} ) n = (t {hH=h} ) p n + = (S h ) p n - (S^Z ~ A 7 nh- 13 / 2 e W {-p 2 nh- 4 ' 3 } . (12) 

This is obtained by taking the h-derivative of the second expression in Prediction 10. 

The asymptotics of (e~ uh ^ u} ' ) )n + = P™ ax (u) as given in Prediction 2 are then obtained by mul¬ 
tiplying (12) by e~ uh , summing over h, and using (18). 


5 Numerical analysis 

We now wish to empirically test the validity of Prediction 2. That is, we wish to investigate 
the behavior of C™ ax (u) and B n (u ) for a variety of values of u, by generating and analyzing the 
sequences up to as large a value of n as possible. 

In full generality then, we will be testing the hypothesis that two one-parameterized sequences 
of positive numbers r n (u), s n (u ), depending on u, have an asymptotic form 

r n (u) ~ A u exp |—An 4 / 7 n 3//7 | n 9 , 
s n (u) ~ A' u exp j— Au 4 ^ 7 n 3,/7 j n 9 , 

respectively, where \, A u , A' u , g, g' are unknown and A ,g,g' are independent of u. (Naturally, we 
will wish to use r n (u ) = C™ ax (w) and s n (u) = B n (u ).) This can also be written as 

log r n (u) = log A u - Au 4 / 7 n 3 / 7 + g logn + o(l). 
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We will make the stronger assumptions that the o(l) is actually 0(n 4 / 7 ) and that a derivative 
form holds, 

/ fiA 7/ 4 / 7 9n\ 

log r n+ 2 (u) - log r n {u) = ( y ^ 4/7 + — J [1 + 0(n -1 )], 

rn ^ = const ■ n 9 ~ 9 + 0(n -11 ^ 7 ). 
s n {u) 

One starts by estimating A using the values of log r n (u) for different values of u and finding the 
slope. 


5.1 Generation of data 

Let b 71t h be the number of bridges of length n spanning a strip of width h. The generating function 
is 

B(z, u) = J2 Khz n e~ uh = Y, B n (u)z n e^ n , 

n,h n 

where B n (u ) is as defined in (3). We have generated all coefficients in B with n < 86 on the square 
lattice. Thus h < 86 also. 

If is the number of upper half-plane SAWs of length n, originating at the origin, with 
maximum height h , the corresponding generating function is 


C™ x (z,u) = Y 


max n —uh 
c n,h z e 


n,h 


Ycr x (n)z n e^, 


where C™ ax (u) is as defined in (2). We have generated all coefficients in (7 max with n < 60 on the 
square lattice. 

The algorithm we use to enumerate SAWs and bridges on the square lattice builds on the 
pioneering work of Enting [8] on self-avoiding polygons extended to walks by Conway, Enting and 
Guttmann [6] with further enhancements over the years by Jensen and others [16, 5, 17]. Below 
we shall only briefly outline the basics of the algorithm and describe the changes made for the 
particular problem studied in this work. 

The first terms in the series for the SAW generating function can be calculated using transfer 
matrix techniques to count the number of SAWs in rectangles of width w and length l vertices long. 
Any SAW spanning such a rectangle has length at least w + l — 2. By adding the contributions from 
all rectangles of width w < n + 1 and length w < l < n — w + l the number of SAW is obtained 
correctly up to length n. 

The generating function for rectangles with fixed width h are calculated using transfer matrix 
(TM) techniques. The most efficient implementation of the TM algorithm generally involves cutting 
the finite lattice with a line and moving this cut-line in such a way as to build up the lattice vertex 
by vertex. Formally a SAW can be viewed as a sub-graph of the square lattice such that every 
vertex of the sub-graph has degree 0 or 2 (the vertex is empty or the SAW passes through the 
vertex) excepts for two vertices of degree 1 (the start- or end-point of the SAW). Thus if we draw 
a SAW and then cut it by a line we observe that the partial SAW to the left and right of this 
line consists of a number of arcs connecting two edges on the intersection, and pieces which are 
connected to only one edge (we call these free ends). The other end of a free piece is either the 
start- or end-point of the SAW so there are at most two free ends. 
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Each end of an arc is assigned one of two labels depending on whether it is the lower end or 
the upper end of an arc. Each configuration along the cut-line can thus be represented by a set of 
edge states {< 7 j }, where 


0 empty edge, 

1 lower arc-end, - . 

Gl — r) J (13) 

2 upper arc-end. 

3 free end. 

Since crossings aren’t permitted this encoding uniquely describes which loop ends are connected. 

The sum over all contributing graphs is calculated as the cut-line is moved through the lattice. 
For each configuration of occupied or empty edges along the intersection we maintain a generating 
function Gs(z) for partial walks with signature S, where Gs(z) is a polynomial truncated at degree 
n. In a TM update each source signature S (before the cut-line is moved) gives rise to a few new 

target signatures S' (after the move of the boundary line) and n = 0,1 or 2 new edges are inserted 

leading to the update Gs’(z) = Gs'(z) + z n Gs(z). Once a signature S has been processed it can 
be discarded. 

Some changes to the algorithm described in [16, 17] are required in order to enumerate the 
restricted SAWs and bridges. We used the recently developed version of the TM algorithm [5, 17] 
in which edge states describe how the set of occupied edges along the cut-line are to be connected 
to the right of the cut-line, i.e, how edges must connect as the cut-line is moved in the transfer 
direction from left to right. 

5.1.1 Further details for SAWs 

Grafting the SAW to the wall can be achieved by forcing one of the free ends (the start-point) to lie 
on the bottom side of the rectangle. In enumerations of unrestricted SAWs one can use symmetry 
to restrict the TM calculations to rectangles with w < n/2 + 1 and l > w by counting contributions 
for rectangles with l > w twice. The grafting of the start-point to the wall breaks the symmetry 
and we have to consider all rectangles with w < n + 1. The number of configurations one need 
consider grows exponentially with w. Hence one wants to minimise the length of the cut-line. To 
achieve this the TM calculation on the set of rectangles is broken into two sub-sets with l > w 
and l < w, respectively. In the calculations for the sub-set with l < w the cut-line is chosen to 
be horizontal (rather than vertical) so it cuts across at most l + 1 edges. Alternatively, one may 
view the calculation for the second sub-set as a TM algorithm for SAWs with its start-point on the 
left-most border of the rectangle. 

In order to measure the maximum height of the SAW all we need to do is extract this information 
from the finite-lattice data set. From the calculation with l > w, where one end-point must lie 
on the bottom of the rectangle, the maximum height is simply w. Similarly for l < w, where 
one end-point must lie on the left of the rectangle, the maximum height is l. The final series is 
then obtained by combining the results from the two cases summing over all possible sizes of the 
rectangles. 

5.1.2 Further details for bridges 

As for SAWs the calculation for bridges is broken into two sub-sets corresponding to bridges which 
span the rectangle from bottom-to-top when l > w = h or from left-to-right when h = l < w. 
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In this case there are some further restrictions on the permissible types of configurations which 
makes the algorithm a little more efficient. The most important is that at all stages during the TM 
calculation a free end must have access to the boundaries of the rectangle otherwise the partially 
constructed SAW would be unable to span the rectangle and hence could not lead to a bridge. So 
any configuration where a free end is embedded inside arcs is forbidden. 

Clearly concatenating two bridges of height i and j gives a bridge of height i + j (we place the 
origin of the second walk on top of the end-point of the first walk). This means that any bridge 
can be decomposed into irreducible bridges, i.e., bridges which cannot be decomposed further, and 
we use a n to denote the number of n-step irreducible bridges. So the generating function B(z,u ) 
for bridges is simply related to the generating function for irreducible bridges A(z, u) 


B(z , u) 


A(z, u) 

1 — A(z, u ) 


(14) 


This fact can be used to extend the series for bridges by the following simple observation. If we 
calculate the number of bridges in rectangles up to some maximal width w max the series for bridges 
will be correct to order 2 w max + 1 but if we do the same for irreducible bridges the series will be 
correct to order 3u> max + 2. Thankfully the number of irreducible bridges can easily be obtained 
from the number of bridges. Consider the number of bridges b n ^ and irreducible bridges a n ^ of 
length n and height h with associated generating functions B k (z) and A k (z). Since a bridge is 
either irreducible or the concatenation of a bridge with an irreducible bridge we get 


and thus 


h-l 

B h (z) = A h (z) + ^2 A h - k {z)B k (z) 
k =1 


h-l 

Ah(z) = B h (z) - YA h _ k {z)B k {z), 

k= i 

which allows us to obtain all generating functions A k (z) recursively from Bh(z) for 1 < h < w max 
and once these are known we can extend the series B(z, y ) by (14) to order n = 3 w max + 2 for any 
1 < h < n. 

We calculated series for Bh(z) up to h = 28 and thus n = 86 which on its own would give a 
series for bridges to order 57 but extracting the data for irreducible bridges and using (14) the 
series to order 86 was obtained. The calculation used a total of some 25,000 CPU hours. The most 
demanding part was the calculation for w = 28 where we used 640 processor cores with 3GB of 
memory used per core. 

5.2 Analysis of data 

The sequences we now analyze are the coefficients of B(z,u) and C max (z, u); that is, {B n (u)e^ n } n 
and {C'™ ax (ii)e^ n } n . For brevity, in the remainder of this paper we will use b n = B n {u)e^ n and 
c+ = C'™ ax (u)e /3n , and will use a n as a placeholder sequence which could stand for either. 

As discussed in [10], if the asymptotic form of a sequence {a n } is 

a n ~ const ■ e^ n ■ Y ■ n 7_1 , (15) 
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2.5 - 2.5 





Figure 1: Bridge series ratios, u = log 2. The top left-hand corner has ordinate e 


then the ratio of successive coefficients is 


r n = 


= 


0"n— 1 


, cr log ii\ g a 2 log 2 iii „ 

1 + —+ - + —+ o 

n 1 a n n. 


2-o 


n 


(16) 


So a ratio plot against n a ~ l should be linear, with perhaps some low-n curvature induced by the 
presence of the 0(1/n) term. 

To illustrate this, we take the generating function for bridges with u = log 2 as a representative 
value, not subject to cross-over effects present near u = 0 and u = oo. Plotting the ratios of 
pushed bridges against 1/n, as shown in Figure la gives a plot displaying considerable curvature. 
Note that the limiting value of as n —> oo is given by the top left corner of the plot. Plotting 
the ratios against l/ra 2 / 3 , as shown in Figure lb gives a plot with reduced, but still significant, 
curvature. In Figures lc and Id we show ratio plots against 1/y/n and against 1/n 4 /" respectively. 
Both are visually linear, but more careful extrapolation reveals that plot lc extrapolates, as n —> 0, 
to a value slightly greater than the numerical value of e^, while plot Id extrapolates to a value 
indistinguishable from the numerical value of e^. Thus this simple sequence of ratio plots gives 
good evidence that the ratios are linear when plotted against 1/n 9 , where 0.5 < 9 < 0.6, and is 
totally consistent with the heuristic expectation, 9 = 4/7 = 0.57142... 

The situation is similar for pushed SAWs (figures not shown), except that for SAWs the possible 
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(a) Plot of estimates of 1 — a against 4. (b) Plot of estimates of a against 4. 

Figure 2: Two different estimators of exponent a. Both are consistent with a = 3/7 = 0.42857.... 


range of values for 6, as chosen by visual linearity of the ratio plots, is even wider, and we estimate 
0.6 < 9 < 0.3, again consistent with the same expected value, 6 = 4/7 = 0.57142... 

In order to more directly estimate the value of the exponent a, we note from (16) that 

[r n ■ e _/3 — 1) ~ const ■ n CT_1 . 

The plot of log(l — r n ■ e ~' 3 ) against logn should be linear, with gradient a — 1. To the naked eye 
it is, so we don’t show it. However there is a small degree of curvature, so we extrapolate the local 
gradients, defined as 

log (l - r n ■ e~P) - log (l - r n -i • e _/3 ) 

1 &n — . . . . , (17) 

log (n — 1) — log(n) 

against 1/n. The results are shown in Figure 2a. The ordinates are estimators of 1 — a, and it is 
clear that they are consistent with the conjectured u = 3/7. 

A second estimator of o can be constructed from (15) by noting that 


log 


log (a n e /3n ) 

= ulogn 

[i + o( log „”)l 

V / 


\ n a J 


As with the previous estimator, the log-log plot of log (a n e -/3n ) against n should be linear, with 
gradient o. We perform this calculation using bridges with u = log 2. To the naked eye the plot is 
linear, so we don’t show it. As in the previous case, there is a small degree of curvature, so we again 
extrapolate the local ratios, as defined above mutatis mutandis, and the resulting extrapolants are 
shown plotted against 1/n in Figure 2b. This is also entirely consistent with the expected value 
3/7. 

We have repeated this calculation for bridges with other values of u, and find that this exponent 
estimate is numerically robust for u in the range e~ u E [0.2,0.8]. Outside this range we see crossover 
effects to other behaviour, as appropriate for u = 0 and u = oo. We have also repeated the 
calculation for SAWs with u = log 2 (and other values of u ), and find good numerical evidence for 
(j « 0.43, supporting the expectation that a = 3/7 exactly for both models. 
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5.3 Estimation of Hi- 

We now take as known the value of the exponent a, and, to very high precision, the value of the 
leading growth constant . We next estimate the value of the (u-dependent) subdominant growth 
constant Hi- Recall that, for random walks (6), the (/-dependence is given by 

log^i = — \ u = const ■ u l - a . 


(The same also occurs for compressed Dyck paths [10].) We find compelling evidence for similar u- 
dependence in the case of both pushed bridges and pushed SAWs. We also find compelling evidence 
that, unsurprisingly, /Ui(SAWs) = Hi (bridges). 

We used several methods to estimate the value of hi. In all cases we utilised both the assumed 
value of the growth constant and the exponent o = 3/7. No one method was especially suc¬ 
cessful for all values of u. As is common with series analysis, all methods involved some sort of 
extrapolation. It was frequently difficult to form a view as to the limiting extrapolant in some 
cases, usually due to a sequence of estimates having a turning point near the end of the range. 
Nevertheless, by using several methods, and looking at the spread between estimates, we were able 
to make moderately precise estimates. 

The methods we used were the following: (i) We first formed the logarithm of the ratio of 
successive coefficients (or, equivalently, the logarithm of the square root of the ratio of successive 
coefficients, ^/\a n /a n - 2), which reduces odd-even oscillations in the ratios). 


log 


CLr 


tin— 1 


— log e p = 


p _ o'log /R , 9 , o- 2 log- m 


n 


1 — <7 


H-h 

n 


2 n 2 ~ 2a 


+ 0(n- 2 - ff ), 


in the case of both sequences a n = b n and a n = c+. We calculated the series expansion of the left- 
hand side of the above equation, known to order n max = 86 in the case of bridges, and n max = 60 
in the case of SAWs. We fitted the terms from order n m i n to order n max to the first two, and then 
the first three terms in the asymptotic expansion, using the Maple procedure NonlinearFit. In 
each case we steadily increased the value of n m i n from (n max — 30) to (n max — 4). (Attempts using 
the first four terms were unsuccessful, due to lack of convergence). 

(ii) The second method involved directly fitting to the coefficients. Since 


log a n — n ■ (3 = n a log Hi +9 log n + const. + O (n <T ) , 


we could fit successive coefficient triples (ak-i, cr Ufc+i), with k increasing up to n max — 1, esti¬ 
mating the three unknowns, log/Ui, g, const, from the three coefficients. We also fitted successive 
coefficient pairs (afc-i, a&) to estimate the first two unknowns. 

The third method (iii) used was a blend of the previous two methods, in that we fitted the 
logarithm of the ratio of successive coefficients to its asymptotic form by the second method - that 
is, a fit using successive pairs or triple of coefficients. 

The next two methods, (iv) and (v) involved estimating hi directly from the ratio of successive 
coefficients, as given by (16). Since 


Sn — Tr 




-1 = 


O' log Ml g 0-- log- Hi 
n n 


n 


1 — (7 


2-2(7 


+ 0 


2-0- 


n 


a plot of s n against n a 1 should be linear, with slope o log h 1 • Slight curvature was evident, so we 
extrapolated the local gradients. 
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In an attempt to get more rapid convergence, we next, (v), eliminated the term 0(1/n) by 
forming the sequence 


t-n — ^ s n (n 1) • s n — i 


0^1 + <j 2 {2<j - i) w m + o (_LA), 

n 1 ~ a n 2-2a y n 2 -a J J 


As above, a plot of t n against n a ~ l should be linear, with slope cx 2 log/ii. Again, slight curvature 
was evident, so we again extrapolated the local gradients. 

The final method, (vi), we used was based on the observation that 


d n = B n (u) n " = fii 1 + 


9 log n 


n u 


+ o{ 


n 


’>) 


so d n provides a sequence of estimators of g\ which can be extrapolated. 

We show in Table 1 the estimates of log /q obtained by the six methods described. The actual 
analysis comprises dozens of pages, so we only give a summary. Rather than attempt to give 
precise error bars, which in any series analysis are perhaps better described as confidence limits, 
we consider it more meaningful to quote parameters to a level of precision such that we expect 
uncertainties to be confined to the last quoted digit. From the various estimates, we first ignore 
obvious outliers, then simply take the average of the remaining estimates, and give this as our 
combined estimate. Comparing results for SAWs and bridges, it appears that g\ is the same, up 
to uncertainty in the quoted best estimates. 

A little numerical experimentation shows that the n-dependence is indistinguishable from log oc 
•u 4 /'. We note that, like random walks (6) and Dyck paths [10], this matches log/i oc v}~ a . Fitting 
the constant of proportionality to the combined estimates leads us to the result 

log g i « —2.62ri 4 / 7 . 


Experimentation with data sets for other u -values is consistent with this result. In order to compare 
with the estimates given below, this formula gives log/ri = —2.913, —2.125 — 1.453 for e~ u = 
0.3, 0.5, 0.7 respectively. 


5.4 Estimation of the exponent g. 

In the previous section we estimated the value of g\ from the series assuming the values of and a. 
Attempting to do the same in order to estimate the exponent characterising the sub-sub-dominant 
term n 9 making the same assumptions was unsuccessful. That is to say, we could not come up with 
any method of analysis that gave a consistent estimate of g for the SAW or bridge case. Note that 
we expect g to be ^-independent, but not necessarily to be the same for SAWs and bridges. 

However if we accept the conjecture (above) that /q (SAWs) = g\ (bridges), one can estimate 
QSAWs — Abridges by studying the asymptotics of the series formed by taking the term by term 
quotient of the SAW and the bridge series. Defining v n = c+/6 n , we expect from the asymptotics 
that v n ~ const ■ n", where a = gsAWs ~ Abridges- 

For e~ u in the range [0.1, 0.6] we find estimates of a, obtained by calculating the gradient of a 
simple ratio plot of r v (n) = v n /v n -i against 1/n, to lie consistently in the range 0.61 to 0.63. For 
e~ u > 0.7 the exponent estimators clearly have a turning point at a value beyond the number of 
coefficients we have available. Accordingly, we cannot estimate a in this tt-regime. 
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Method 

CO 

o 

II 

1 

e~ u = 0.5 

e~ u = 0.7 

(i) Bridges, 2 terms 

> -2.98 

ps -2.177 

> -1.51 

(i) Bridges, 3 terms 

ps -2.89 

> -2.180 

< -1.30 

(ii) Bridges, 2 terms 

-3.0 

-2.15 

-1.5 

(ii) Bridges, 3 terms 

-2.92 

n.c. 

-1.46 

(iii) Bridges, 2 terms 

-2.905 

-2.13 

-1.47 

(iii) Bridges, 3 terms 

n.c. 

n.c. 

-1.46 

(iv) Bridges 

-2.96 

-2.14 

-1.47 

(v) Bridges 

-2.92 

-2.14 

n.c. 

(vi) Bridges 

-3.0 

-2.21 

-1.5 

Combined Estimate 

-2.92 

-2.14 

-1.465 

(i) SAWs, 2 terms 

> -3.00 

> -2.21 

> -1.55 

(i) SAWs, 3 terms 

n.c. 

ps -2.14 

> -1.485 

(ii) SAWs, 2 terms 

> -2.98 

-2.13 

-1.5 

(ii) SAWs, 3 terms 

< -2.90. 

n.c. 

-1.45 

(iii) SAWs, 2 terms 

-2.87 

-2.135 

-1.47 

(iii) SAWs, 3 terms 

-2.92. 

n.c. 

-1.45 

(iv) SAWs 

-2.84 

-2.12 

-1.42 

(v) SAWs 

-2.89 

-2.15 

-1.47 

(vi) SAWs 

n.c. 

n.c 

n.c. 

Combined Estimate 

-2.91 

-2.13 

-1.45 


Table 1: Estimates of log/Ui for SAWs and bridges for various e u values, using the six methods 
described. ( n.c. means no convergence). 
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In order to estimate the individual exponents, we re-analyse the original series using not only 
our best estimate of and the value a = 3/7, but also the estimate of fi\ obtained in the previous 
section. We can, using these values, divide the original series coefficients by e n ^ ■ /ri(u) n<T . The 
resultant coefficients should behave asymptotically as const ■ n 9 . Thus g can be estimated in a 
variety of ways, but given the uncertainty in the value of ji \, there is no point in using sophisticated 
methods. Rather, we just look at the gradient of the ratio plot of the divided coefficents, when 
plotted against 1/n. The ratio of succesive coefficients in this case should behave as l+g/n+o(g/n). 
Thus plotting the ratios against 1/n should give a straight line with gradient g. 

In this way we estimate (/saws ~ 0.15, 0.15, 0.16 for e~ u = 0.3, 0.5, 0.7 respectively, and 
Abridges ~ —0.5, —0.5, —0.44 for e~ u = 0.3, 0.5, 0.7 respectively. These results are reasonably 
consistent with the preceding estimate that the gap between the exponents should be in the range 
0.61 to 0.63. Our estimate for (/bridges has an uncertainty of about ±0.05, while that of gsAWs is 
±0.02, so our estimates have a difference of 0.65 ± 0.07, which is reasonably consistent with the 
direct estimate found above. 

The heuristic arguments in the previous section predict that (/bridges = — = —0.4642857..., 

in agreement with our analysis above, while (/saws = fg = 0.18750, again in agreement with the 
above analysis. Similarly, a = ^ = 0.651785 ..., as against our direct estimates of a « 0.65 and 
a ~ 0.62. 

6 Conclusion 

We have studied several models of self-avoiding walks (SAWs) and polygons (SAPs) which are 
constrained to lie in the upper half-plane and are subjected to a compressive force. The force is 
applied to the vertex or vertices of the walk located at the maximum distance above the boundary 
of the half-space. We have in particular focused on three types of objects: SAWs, SAPs and self¬ 
avoiding bridges. In each case, we have considered the partition function of objects of size n, with 
the aim of determining the asymptotic behaviour of these partition functions in the n —> oo limit. 

We used the conjectured relation with the Schramm-Loewner evolution to predict the asymp¬ 
totic forms of the partition functions, including the values of the exponents. These values, as stated 
in Prediction 2, are that (for SAWs, bridges and SAPs respectively) 

C™ x (u) ~ A+ n 3/16 exp{-Am 4/7 n 3/7 }, 

B n (u) ~ A b u n~ 13 / 28 exp{—A lU i/7 n 3/7 }, 

P™ x {u) ~ AP+n- 11 ^ 7 exp{-A 2 u 4 / 7 n 3 / 7 }, 

where u > 0 corresponds to the compressive force regime, Ai, A 2 > 0 are constants and A+, A^, Au + 
are functions of u. 

Finally, we tested the predictions for SAWs and bridges by analysing exact enumeration data 
and found them to be in agreement (within the range of uncertainty resulting from analysis of 
limited series data). A similar analysis for SAPs has been performed elsewhere. 
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A Asymptotics of an integral 

Several times in this paper, we used the asymptotics of a particular integral. Here we justify that 
relation. 

Lemma 1. If a, 6, k > 0, r G M, then as m —> oo, 

km 


f ( km \ ) 2r+l r 1 'l 

J y r exp I — I ^ + byJ > dy ~ C m 2 ( 1+a ) exp Am 1 + a j , (18) 


where 


c = 6~( 1+r ) 27r 


1 + a 


A = 


(a + l)^^^ (a + 1) f k \ !+ 


a 


a 


— 6 !+° 
6 ° 


Proof. Using the change of variables x = by, we see it suffices to prove the result for 6=1. We 
write km = n l+a , so it suffices to show that 


x exp 


1+a 


n 


ax'- 


+ x ) > dx = 


2 vr 

1 H - CX 




i + o^-v 2 ) 


Using the substitution x = n l ^ 2 y + n, 

,1+a 


x r exp 


n 




+ x > dx 


= n r+ 3 


*jC.( 


1 H—exp < —n 
n 


n 


+ ~ 1 = + 1 ) } dy 


= n r+ 2 e _n ( 1+ «) 


a(n 1 / 2 y + ro)" 

1 + ( y/Vn)) r exp {-n#(y/Vn) } ciy, 


where 

s(*) = 5a(*) = « _1 ((i + 1)"“ - l) + 1. 
Using the substitution t = y/\/n, 

,1+a 


n -(r+i) e n(l+i) 


x r exp 


n 


as" 


+ x > dx = \m / (1 + t) r exp {—ng(f)} dt. 


'-1 


Since g(0) = 0,g'(0) = 0,g"(0) = a + 1, we have 

g{t) = + o(f), 

which is valid for |t| < 1/2. Note that there exists c > 0 such that 

g(t) > cn ~ 2 / 3 , |t| > n" 1 / 3 , 
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t > 0 


g(t) > t — a 1 _ , 

1 Oi T 1 1 

g(t)> - ---, 

yw ~ (1 + t)° a ~ 2 


Using these estimates in order, we see that 


j | ^l + t) r exp{-ng(t)}l {W > n -i /3} dt = o(n\ rl exp|-cn 1/3 |) = o(n 1/2 ), 

poo poo 

/ (1 + t) r exp{— ng(t)} dt < e n/ " / (1 + t) r exp{— nt} dt < 0(n^ e~ n ^ a ) = o(n -1 / 2 ), 

J2/a J2/a 

r—l +n -1 

J (1 + t) r exp{— ng(t)} dt = o(n -1 / 2 ). 


Hence, 


/ OO /■)! 1//3 

(1 + t) r exp {—rag'(t)} (it = o(n -1 ^ 2 ) + \/n / (l + i) r exp {—n g(i)} dt. 

-l 7-n - 1 / 3 

For |i| < n _1//3 , we write 

(1 + ty e ~ n = [1 + O(t)] exp j-n i 2 + 0(t 3 )) | 

= e - n (“ +1 ) t2 / 2 [1 + O(t) + 0(ni 3 )]. 

Here we have restricted to |i| < rU 1 / 3 so that nt 5 = 0(1) and hence the expansion e°( n * 3 ) = 
1 + 0(nt 3 ) is valid. 

Then 


/•n -1 / 3 .n 1 / 6 

/ e ~n(a+l)t 2 /2 dt = / e _( Q+ l )s 2/2 = 

' — n~ V 3 «/—n 1 / 6 


(o; + 1 ) 


-o(n 1/2 ), 


/»n - 1 / 3 

n (\t\+n\t\ 3 ) e~ n ^ a+1)t2/2 dt< 

J-n- 1 / 3 


Isl + Isr 


e ~(a+l)s 2 /2 — 0(n~ 1//2 ). 
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